arXiv:gr-qc/9706010vl 3 Jun 1997 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 February 2008 


Pulsation modes for increasingly relativistic polytropes 

Nils Andersson^ and Kostas D. Kokkotas^ 

^ Department of Physics, Washington University, St Louis MO 63130, USA 
^ Department of Physics, Aristotle University of Thessaloniki, Thessaloniki SfOOO, Greece 


Accepted 1997 (?). Received 1997 (?); in original form 1997 


ABSTRACT 

We present the results of a numerical study of the fluid /, p and the gravitational 
w modes for increasingly relativistic nonrotating polytropes. The results for / 
and tc-modes are in good agreement with previous data for uniform density 
stars, which supports an understanding of the nature of the gravitational wave 
modes based on the uniform density data. We show that the p-modes can 
become extremely long-lived for some relativistic stars. This effect is attributed 
to the change in the perturbed density distribution as the star becomes more 
compact. 
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1 INTRODUCTION 


In a recent paper (|Andersson, Kojima fc Kokkotas 1996|) we presented a detailed survey 
of the pulsation modes of the simplest conceivable stellar model: a nonrotating star with 
uniform density. The reasons for chosing this, admittedly unrealistic, model were twofold: 
First of all, the analytic solution to the TOV equations for uniform density is well known 
( ^chutz 1985|) . This means that calculations could readily be done for many different values 
stellar compactness and average density. Secondly, more realistic models for a neutron star 
are “almost constant density” so the results for this simple model may not be too different 
from what one would hnd for realistic equations of state. 

The main conclusions of our study were: 


• that gravitational-wave {w) modes ([Kokkotas and Schutz 1992|) exist for all stellar 
models, even though the modes become very short lived for less relativistic stellar models. 









2 Nils Andersson and Kostas D. Kokkotas 


Moreover, the w-mode spectra are qualitatively similar for axial and polar perturbations 
(for a description of the two classes of relativistic perturbations - also referred to as odd and 
even parity perturbations - see Thorne and Campolattaro (|1967| ) or Kojima (|1992| )). Since 
axial perturbations do not couple to pulsations in the stellar fluid ( [hhorne fc Campolattaro 
19671) we concluded that the tc-modes are pure “spacetime” modes that do not depend on 


the particulars of the fluid for their existence. This notion is supported by results obtained 
for various model scenarios ([Andersson, Kokkotas fc Schutz 1996| ; [Andersson 1996|) . 

• that the behaviour of the various pulsation modes changes dramatically as the star 
is made more compact than R ^ 3M (we use geometrized units c = G = 1). For such 
ultracompact stars the existence of a peak in the effective curvature potential outside the 
surface of the star will affect the various modes. Some modes can basically be considered as 
“trapped in a potential well”, cf. Chandrasekhar and Ferrari ( p.991[) and Kokkotas ([1994[) . 
This trapping has the effect that some w-modes become extremely long-lived. 

• that there are “avoided crossings” between the fluid /-mode and the various polar w- 
modes for extremely compact stars. This suggests that the /-mode should be considered as 
the hrst in the sequence of trapped modes for these stars. 


The results of the mode-survey for uniform density models helped improve our under¬ 
standing of the origin of, and the relation between, various pulsation modes of relativistic 
stars. We would expect the main conclusions to hold also for more realistic equations of state 
(at least qualitatively), but the motivation for conhrming these expectations by actual cal¬ 
culations is strong. For example, axial tc-modes have so far only been calculated for uniform 
density stars. Although our present understanding of the origin of the tc-modes indicates 
that axial modes must exist for ^ stellar models it is important to verify that this is, in¬ 
deed, the case. Furthermore, more realistic stellar models support other families of modes. 
Of particular interest for gravitational-wave physics are the p-modes (that correspond to 
pressure waves in the stellar fluid). The fact that the properties of the /- and tc-modes 
change dramatically as the star becomes increasingly relativistic leads to questions whether 
the p-modes are also affected in an interesting way. With this short paper we address these 
issues. 

We describe results obtained for the polytropic equation of state 

p = Kp^ , 

where k = 100 km^. Our main study was for F = 2, but we also considered other (especially 
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larger) values of F. Our results are basically an extension of those presented by Andersson, 
Kokkotas and Schutz (|1995| ). A detailed description of the way that we extract the complex 
frequencies of the various pulsation modes can be found in that paper. In order to keep the 
present paper short, we will discuss neither this numerical approach nor the various pertur¬ 
bation equations here. The equations governing polar perturbations are described in detail 


by Lindblom and Detweiler (|1981| ), while the axial equations are given by Chandrasekhar 


and Ferrari (|1991|) . 


2 NUMERICAL RESULTS 

The results of our investigation for increasingly relativistic polytropes can be summarized 
in a single figure. In Figure |I| we show the inverse damping rate (represented by Im uM) of 
each mode as a function of the pulsation frequency (Re utM) for a F = 2 polytrope. As the 
central density of the stellar model is varied each mode traces out a curve in the complex 
(u-plane. In the hgure we have indicated (by diamonds) the densest polytropic model that is 
stable to radial perturbations. That is, the specific model for which the mass (M) reaches a 
maximum as a function of the central density {pc)- For the F = 2 polytrope used to calculate 
the data in Figure |^the marginally stable model has central density pc = 5.7 x 10^® g/cm . 

It is clear that Figure |I] contains a considerable amount of information. We now proceed 
to discuss the results for each separate class of modes in more detail. 


2.1 Gravitational wave modes 


As far as the tc-modes are concerned our investigation does not add much to the study for 
uniform density stars ( Andersson, Kojima fc Kokkotas 1996 ). The results for polytropes 
are, both qualitatively and quantitatively, in excellent agreement with those for the uniform 
density model. From the upper panel in Figure it is clear that both the axial and the polar 
tc-modes can be divided into two separate families. The two families are distinguished by 
the behaviour of the mode-frequencies as the star becomes very relativistic. For one family 
of modes - the interface w-modes ( [Andersson, Kojima fc Kokkotas 1996^ [Leins, Nollertj 
and Sohel 1993[) - the frequencies approach constant (complex) values for large pc- For the 
specific polytrope under consideration we find that the mode-frequencies hardly change at 
all for pc > 1.9 X 10^® g/cvsi . The modes in the second family - the curvature tc-modes 
- are characterized by an increase in |a;M| with pc for less compact stars, but before the 
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star becomes radially unstable \u}M\ reaches a maximum value. The modes in this family 
then become extremely slowly damped as pc increases further. The slowest damped w-modes 
become long-lived as R < 3.5M or so. The decreased damping rate is most likely due to 
the increased influence of the curvature potential barrier (the peak of which is at i? ~ 3M). 
However, it is important to realize that the T = 2 polytrope that we consider is always less 
compact than R = 3.2M. Hence, the surface of these stars is never located inside the peak 
of the curvature potential. That the modes still become extremely slowly damped shows 
that the exterior potential can “trap” gravitational waves effectively even though there is 
no “potential well inside a barrier” for these stars, cf. the discussion by Detweiler (|1975|) . 

As already mentioned, the present results are very similar to those for uniform density 
stars. Especially worth noticing are the similarities between the axial and the polar tc-modes 
in Figure 0. Since the axial modes cannot induce oscillations in the stellar fluid, this indicates 
that the character of the tc-mode depends solely on the properties of the curved spacetime. 
In other words: the tc-modes are “spacetime” modes the details of which do not depend on 
the dynamics of the stellar fluid (|Andersson, Kojima fc Kokkotas 1996| ; [Andersson, Kokkotas 
fc Schutz 1996| ). 


2.2 Fluid modes 

The present results for the fluid /-mode are also in good agreement with the uniform density 
study ( Andersson, Kojima fc Kokkotas 1996| ). For less compact stars the damping rate of 
the /-mode increases with p^. This simply means that the star radiates gravitational waves 
more efficiently as it becomes more compact. But, in a similar way to the curvature w- 
modes, the /-mode becomes slower damped once the star approaches R = 3M. Again, this 
can be understood as the increasing influence of the curvature potential (Petweiler 1975] ; 
Andersson, Kojima fc Kokkotas 1996|) . 

A notable difference between the results for the T = 2 polytrope and the uniform density 
star is the absence of “avoided crossings” between the polar w-modes and the /-mode, cf. 
Figure 2 of Andersson, Kojima and Kokkotas (|1996|). The reason for this is that the T = 2 


polytrope never becomes sufficiently compact for such avoided crossings to occur. For the 
uniform density model the hrst avoided crossing occured for R ^ 2.3M, and for the T = 2 
polytrope we always have R > 3.2M. One would, however, expect avoided mode-crossings 
to exist for polytropes that can be made very compact. We have verihed the existense of 
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avoided crossings for a F = 5 model (which can be made as compact as i? ~ 2.39M). 
The existence of avoided crossings thus seems to be a generic feature of extremely compact 
relativistic stars. 

The results we present for p-modes are both new and remarkable. Fluid p-modes have not 
previously been calculated for extremely relativistic stellar models. As can be seen in Figure]^ 
these modes behave in a peculiar ways as Pc increases. While the behaviour for less relativistic 
stars is analogous to the /- and the w-modes - the damping rate increases when the star 
becomes more compact - it is different for the very relativistic models. Once the star has 
become radially unstable |a;M| attains a maximum. Then the pulsation frequency Re uM 
typically decreases monotonically. At the same time the damping rate decreases drastically 
and Im ujM has a sharp minimum. For some modes we hnd several such minima of Im ujM, 
cf. Figure |I[ The corresponding values of the imaginary part of the mode-frequency are 
actually so small that our code does not have sufficient precision to distinguish them from 
zero. It is, of course, important to establish that the sign of Im ojM does not change. A 
change in sign would indicate that the mode becomes linearly unstable. We believe that our 
calculations were sufficently accurate to establish that this does not happen ~ the p-modes 
are all stable. But it is interesting to note that they become extremey long-lived (very poor 
radiators of gravitational waves) for some values of the central density. 

How can we understand the peculiar behaviour of the p-mode frequencies? Let us hrst 
consider the pulsation frequency Re ujM, and compare our results to ones for simpler stellar 
models. The well-known result for the /-mode, established for incompressible fluid spheres 
by Kelvin in 1863 ([Tassoul 1978|) , suggests that our mode should approach 



for less relativistic stars. A similar result, for compressible homogeneous spheres ([Tassoul 


1978|) , shows that the p-modes ought to approach 



p — [An -h R A2 -K /(/ 1)] ( — ) , (2) 


where 


2An = [2/ + 3 + n{2n + 21 + 5)]F - 4 , n = 0,1, 2,... 


( 3 ) 


and F is the polytropic index. As can be seen from Figure our numerical results agree quite 
well with these approximations. That is, the pulsation frequency of each fluid mode changes 
in the anticipated way as the star become more relativistic. Alternatively, the results in 
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Figure ^ can be taken as justification for using the approximate results also for relativistic 
stars. They clearly provide useful estimates also for objects of neutron star compactness, 
R/M ^ 5. 

The peculiar behaviour of the p-mode damping rates can be explained by changes in 
the matter distribution inside the star. According to the standard quadrupole formula, 
gravitational waves due to the fluid motion in the star can be estimated by 

h ~ 5p(r, t) dr =-/ r^5p{r,uj)dr , (4) 

r JO oR r Jo 

where 6p is the Eulerian variation in the density. We have assumed that the perturbation 
is monochromatic with a harmonic time-dependence 6p{r,t) = 6p{r,u})e^^^. In the notation 
of Lindblom and Detweiler ( |1981|) the required density perturbation follows from (for the 
quadrupole) 


5p{r,u) = —N 




P + P 


eN\M + 47rpr=^)iy - (p + p)Ho 


P 


i-r 


kT 


( 5 ) 


where we have used the fact that the acoustic wave speed in a polytrope is 

To get to a useful result, we must also figure out how each complex-frequency mode 
contributes to the relevant physical quantity (such as the density variation). The reason 
is obvious: For complex frequencies a; the 5p(r, cj) that follows from (|]) will be complex, 
whereas the corresponding physical quantity should be real valued. The analysis of this 
problem proceeds exactly as for the analogous black-hole problem ([Andersson 19971) . The 
main steps are i) note that the modes come in pairs and —cj* (where the asterisk de¬ 
notes complex conjugation) ii) Since all perturbation equations [see Lindblom and Detweiler 
(p.981|)] contain only o;^, and (o;^)* = (—o;*)^, we can conclude that 


W’’.-O = ■ 


( 6 ) 


From a specific mode frequency a;„ we therefore get a contribution (for more details see 
Andersson ( j 1997]) ) 


6p{r,t) rsj 2Re \6p{r,uJn)e 


t 


( 7 ) 


Both the real and the imaginary part of the eigenfunction 6p{r,LJn) will thus be important. 

Let us now assume that the gravitational-wave damping is sufficiently slow that the mode 

is essentially undamped for an entire cycle (this is a reasonable assumption for the slowly 

damped fluid modes). Averaging over one period we get 

2 


< e 


—2Ima;nt 


r'^(5p(r, cOnjdr 


= e 
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e{uJn) 


( 8 ) 
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The value of e(a;„) provides a measure of how “efficient” a specihc mode is as radiator 
according to the quadrupole formula. An association between the p-mode minima of Im uM 
and minima in e(un) would indicate that the features seen in Figure are due to changes 
in the perturbed density distribution as Pc increases. In Figure ^ we compare Im ujM to 
e(a;„). This hgure seems to establish the anticipated correlation between minima in the two 
functions for the first few p-modes. 

Having found the likely explanation for the extremely slow damping of various p-modes 
for some values of pc, it is worthwhile to discuss how the eigenfunction Sp{r, ujn) changes as 
we vary pc. Let us hrst consider the /-mode. This mode is distinguished by the fact that 
the eigenfunctions i) have no nodes inside the star and ii) grow monotonically towards the 
surface of the star. The p-modes are different; in general there are n nodes in both the real 
and the imaginary part of 6p{r, ut) for the nth p-mode. But the location of these nodes change 
as we vary pc. We hnd that the distribution of 6p{r, u) changes drastically with varying pc. 
Specifically, for the hrst p-mode and pc < 1 x 10^® g/cm^ the single node is close to the 
surface of the star and the bulk of 6p{r,uj) is located in the outer 50% of the star. But for 
Pc > 3 X 10^® g/cm^ the node has moved close to the centre of the star and the bulk of 
6p{r, uj) is now located inside half the radius of the star. It is understandable that these two, 
very different, conhgurations can be rather different as radiators of gravitational waves. For 
all p-modes the trend is that the bulk of the eigenfunctions move towards the centre of the 
star as we increase pc- In a way, this means that the p-modes for the very relativistic star 
have the same features as the p-modes ( fTassouI 1978|) . 


3 FINAL COMMENTS 


We have presented numerical results for the pulsation modes of nonrotating relativistic 
polytropes. These results agree well, both qualitatively and quantitatively, with previously 
obtained results for uniform density stars ( Andersson, Kojima fc Kokkotas 1996 ). This 
provides further evidence for our understanding of the nature of the various pulsation modes, 
especially the notion that the gravitational wave tc-modes are mainly due to the properties 
of the curved spacetime ( |Andersson, Kokkotas fc Schutz 1996|) . Our present survey also 
considers the p-modes of a relativistic star. We hnd that the polytrope p-mode frequencies 
are close to an approximation based on the pulsation of a compressible homogeneous sphere 
for less relativistic stars. But we also show that the p-modes change in a peculiar way as 
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the star becomes more compact. Specifically, we find that each p-mode can be very slowly 
damped (in fact, almost undamped) for some values of the central density. We have shown 
that this feature can be understood in terms of the change in the distribution of the perturbed 
density with varying central density. Although the result is not of tremendous astrophysical 
importance (since the affected stars are not stable to radial perturbations, anyway) it is 
interesting to note that the effectivity with which a p-mode radiates gravitational waves can 
vary considerably between two stellar models that have almost identical central densities. 
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Figure 1. The change in the complex mode-frequencies as the central density of the stellar model varies (pc increases in the 
direction of the arrows). The data are for a F = 2 polytrope. The gravitational tf;-modes are shown in the upper panel, while the 
much longer lived fluid /- and p- modes are in the the lower panel (The /-mode is the leftmost mode in the lower panel.). Axial 
tc-modes are represented by dashed lines while polar modes are shown as solid lines. The diamonds on each curve indicates the 
densest stellar model that is stable to radial perturbations. 
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Figure 2. The pulsation frequencies of the fluid modes as functions of the stellar compactness KjM. We compare our results 
for /, po, Pi and p 2 (solid lines, from bottom to top) to approximate results for homogeneous stellar models (dashed lines). 
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Figure 3. We compare the damping rate of the fluid pulsation modes (represented by Im loM and solid curves) to the efficieny 
measure e{LOn) (dashed curves). Both quantities are shown as functions of the central density of the stellar model (in units 
of g/cm^) The scale for Im uM is logarithmic and ranges from 10~® to 10“^ while the scale for e{uJn) is arbitrary. The 
data are for a) the /-mode, b-d) the first three p modes. Of interest here is a possible correlation between minima in the two 
functions. Such a correlation explains the minima in the damping rate of the p-modes in terms of changes in the perturbed 
density distribution in the star. (The /-mode is only included for comparison — there should be no minima in e{ujn) for that 
mode.) 




